Measuring the photon distribution by ON/OFF photodectors 
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Reconstruction of photon statistics of optical states provide fundamental information on the nature of any 
optical field and find various relevant applications. Nevertheless, no detector that can reliably discriminate the 
number of incident photons is available. On the other hand the alternative of reconstructing density matrix 
by quantum tomography leads to various technical difficulties that are particular severe in the pulsed regime 
(where mode matching between signal an local oscillator is very challenging). Even if on/off detectors, as usual 
avalanche PhotoDiodes operating in Geiger mode, seem useless as photo counters, recently it was shown how 
reconstruction of photon statistics is possible by considering a variable quantum efficiency. Here we present ex- 
perimental reconstructions of photon number distributions of both continuous-wave and pulsed light beams in a 
scheme based on on/off avalanche photodetection assisted by maximum-likelihood estimation. Reconstructions 
of the distribution for both semiclassical and quantum states of light (as single photon, coherent, pseudofhermal 
and multithermal states) are reported for single-mode as well as for multimode beams. The stability and good 
accuracy obtained in the reconstruction of these states clearly demonstrate the interesting potentialities of this 
simple technique. 

PACS numbers: 42.50.Ar, 42.50.Dv, 03.65.Wj 
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I. INTRODUCTION 



The evaluation of diagonal elements of the density matrix 
for quantum optical states, i.e. of the statistical distribution of 
the number of photons, provides fundamental information on 
the nature of any optical field and finds various relevant appli- 
cations, ranging from studies on foundations of quantum me- 
chanics 1 1] to quantum information \ 2\, and quantum metrol- 
ogy. Despite the importance of photon distribution, photon 
detectors allowing an effective discrimination among differ- 
ent number of incident photons are not yet available. Among 
the possible candidates to number resolving photo-detectors, 
PhotoMultiplier Tubes (PMT's) |3] and hybrid photodetectors 
0] have the drawback of a low quantum efficiency, since the 
detection starts with the emission of an electron from the pho- 
tocathode. On the other hand, solid state detectors with inter- 
nal gain, in which the nature of the primary detection process 
ensures higher efficiency, are still under development. Highly 
efficient thermal photon counters have also been used. How- 
ever, since they operate in cryogenic conditions they are far 
to allow common use EH. Furthermore, their efficiency is 
limited by the optical window for entering the cryostat. 

Quantum tomography provides an alternative method to 
measure photon number distributions jg]. However, the to- 
mography of a state, which has been applied to several quan- 
tum states 1 10], needs the implementation of homodyne de- 
tection, which in turn requires the appropriate mode matching 
of the signal with a suitable local oscillator at a beam splitter. 
This technique is therefore, in general, not of simple imple- 
mentation and, in particular, such mode matching is a very 



challenging task in the case of pulsed optical fields. 

On the other hand, the photodetectors usually employed in 
quantum optics, such as Avalanche PhotoDiodes (APD's) op- 
erating in the Geiger mode |6, 11] (that have relatively large 
quantum efficiencies), are not suited for distinguishing differ- 
ent number of incident photons, since they have the obvious 
drawback that the breakdown current is independent of the 
number of detected photons. The outcome of these APD's 
is either "off (no photons detected) or "on" i.e. a "click", 
indicating the detection of one or more photons. Actually, 
such an outcome can be provided by any photodetector (PMT, 
hybrid photodetector, cryogenic thermal detector) for which 
the charge contained in dark pulses is definitely below that of 
the output current pulses corresponding to the detection of at 
least one photon. Notice that for most high-gain PMT's the 
anodic pulses corresponding to the "no-photons" ("no-click") 
event can be easily discriminated by a threshold from those 
corresponding to the detection of one or more photons. On 
the other hand, , as we will describe in the next paragraph, 
these detectors can be used for reconstructing photon statis- 
tic when measurements at different quantum efficiencies are 
performed. 

In this paper we present in some details (see |8] for a 
shorter summary of these results) experimental reconstruc- 
tions of photon number distributions of both continuous- wave 
and pulsed light beams in a scheme based on on/off avalanche 
photodetection assisted by maximum-likelihood estimation. 
Reconstructions of the distribution for both semiclassical and 
quantum states of light (as single photon, coherent, pseu- 
dothermal and multithermal states) are reported for single- 
mode as well as for multimode beams. 
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II. THEORETICAL METHOD 

The statistics of the "no-click" and "click" events from an 
on/off detector, assuming no dark counts, is given by 



(1) 



and p>o(v) — 1 — Po(v)' where g n = (n\g\n) is the photon 
distribution of the quantum state g and r\ is the quantum effi- 
ciency of the detector, i.e. the probability of a single photon 
to be revealed. At first sight the statistics of an on/off detector 
appears to provide quite a scarce piece of information about 
the state under investigation. However, if the statistics about 
Po{t]) is collected for a suitably large set of efficiency values 
then the information is enough to reconstruct the whole pho- 
ton distribution g n of the signal, upon a suitable truncation at 
n of the Hilbert space. 

The reconstruction of photon distribution through on/off 
detection at different efficiencies has been analyzed 1 121] and 
its statistical reliability investigated in some details Hal . In 
addition, the case of few and small values of -q fl4ll has been 
addressed. However, whilst these theoretical studies found 
an application to realize a multichannel fiber loop detector 
1 1 611 , an experimental implementation of this technique for 
reconstructing photon distribution of a free-propagating field 
is still missing. In view of the relevance of photon distribution 
for applications in quantum information and foundations of 
quantum mechanics, our purpose is to show that a reconstruc- 
tion of the photon distribution by using this technique can be 
effectively realized gathering results obtained from measure- 
ment at different quantum efficiencies. As we will see this 
method leads to excellent results both for free-propagating 
continuous-wave (cw) and pulsed light beams, for both single- 
mode semiclassical and quantum states, as well as for multi- 
mode states. 

The procedure consists in measuring a given signal by 
on/off detection using different values <q y (v — 1, K) of 
the quantum efficiency. The information provided by ex- 
perimental data is contained in the collection of frequencies 
/„ = /o(f?iy) = na v /n v where uq u is the number of "no 
click" events and n v the total number of runs with quantum 
efficiency r\ v . Then we consider expression Q as a statisti- 
cal model for the parameters g n to be solved by maximum- 
likelihood (ML) estimation. Upon defining p v = po(r)v) 
and A vn — (1 — r\y) n we rewrite expression Q as p„ = 



Since the model is linear and the parameters 



to be estimates are positive (LINPOS problem), then the so- 
lution can be well approximated by using the Expectation- 
Maximization algorithm (EM) [17]. By imposing the restric- 
tion J2 n Qn = 1, we obtain the iterative solution 



= Q 



K 



(2) 



where p„[{f?n }] are the probabilities p v , as calculated by us- 
ing the reconstructed distribution {g^n} at the i-th iteration. 
Before going to the experimental implementation we have 
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FIG. 1: Plot of the total error as a function of the iteration num- 
ber for the various state of the field (corresponding to the states ex- 
perimentally investigated, see next Sections). For the thermal dis- 
tribution the value of the total error is plotted at each iteration step 
(value on the x axis), while for the other states the value of e^*' is the 
average over an suitable set of iterations. 



performed several numerical simulations in order to check the 
accuracy and reliability of this method by varying the dif- 
ferent parameters. Since the solution of the ML estimation 
is obtained iteratively, the most important aspect to keep un- 
der control is its convergence. Of course, the degree of con- 
vergence at a given step can be checked evaluating the log- 
likelihood function C = log L, 

K 



L = Y[ P :°»(i- Pu y 



However, a more suitable parameter is given by the total ab- 
solute error at the i-th iteration i.e. 



K 



v=0 



(3) 



Indeed, the total error measures the distance of the probabili- 
ties Pi, [{£>«■*}], as calculated at the i-fh iteration, from the ac- 
tual experimental frequencies and thus, besides convergence, 
it quantifies how the estimated distribution reproduces the ex- 
perimental data. The total distance is a decreasing function of 
the number of iterations. Its stationary value is proportional 
to the accuracy of the experimental frequencies {/„}. For fi- 
nite data sample this value is of order 1 / y/nj for each value 
of r\ v , giving us a rough estimate of nj \fnZ for the total error 

for the reconstructed probability p v [{Qn }]• If me stationary 
value of eW is of this order we have double checked the con- 
vergence of the whole method. Notice that these properties 
are not shared by the log-likelihood: its stationarity certainly 
reveals convergence, but the value depends on the statistics 
to be retrieved and so can not even be estimated by a priori 
considerations. 

In order to check the convergence of the iterations in (0, we 
run simulated experiments to reconstruct some of the states 
subsequently investigated experimentally (see next Sections 
and Figs. I3I7I >. The results are shown in Fig. [0 where the total 
error (ratio its stationary value) is reported as a function of the 
number of iterations. As it is apparent from the plot the total 
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error shows a transient behaviour and then quickly converges 
to its stationary value. The rate of convergence depends on the 
signal under investigation. Our results show that the iterative 
algorithm always converges and the asymptotic value of 
is of the expected order on magnitude. 

An estimate of the confidence interval on the determination 
of the element g n can be given in terms of the variance 

cr n = {U F n )- X '\ (4) 

M being the total number of measurements, and F n the 
Fisher's information |7] 

F n = ^f^)\ (5) 

^ q V \OQnJ 

where 

Vv _ Z^n A vn Q n 

represents the renormalized probabilities of the no-click event 
at quantum efficiency r\ v . 

III. EXPERIMENTAL DATA IN CW REGIME 

In the following we present various different applications 
of this method both to cw and pulsed regime, with the pur- 
pose of demonstrating the potentialities of this technique. For 
what concerns cw regime we have studied reconstruction of 
diagonal element of the density matrix for single photon Fock 
states and a weak coherent one. 
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FIG. 2: Set up for heralded single PDC photon reconstruction. An 
Ar + laser beam pumps a type II BBO crystal generating collinear 
degenerate PDC. After an anti-UV filter, pairs are split on a PBS. 
Observation of a photon on trigger detector D (preceded by a lens, L, 
and an interference filter IF) starts a TAC ramp which is eventually 
closed by a stop signal deriving from the observation of a photon in 
the second detector. 

The single photon states have been generated by produc- 
ing Parametric Down Conversion (PDC) heralded photons. In 
some more detail (see Fig. |2j, pairs of correlated photons have 
been generated by pumping in collinear degenerate geometry 
a 5 mm long /3-Barium-B orate (BBO) crystal, cut for type II 
phase matching (i.e the two photons of the pair have orthogo- 
nal polarizations), with a 0.3 W cw Argon ion laser beam with 
wavelength of 35 1 nm. The heralded single photon scheme is 



based on the specific properties of PDC emission. PDC is a 
quantum effect without classical counterparts and consists of a 
spontaneous decay, inside a non-linear crystal, of one photon 
from a pump beam (usually generated by a laser) into a couple 
of photons conventionally called signal and idler. This decay 
process obeys (phase matching laws) to energy conservation 

UJQ = Wj + LJ S (7) 

and momentum conservation 

k = h + k s (8) 

where uiq, 0Ji,uj s are the frequencies and fcn, fcj, k s the wave 
vectors of pump, idler and signal photon respectively. Fur- 
thermore, the two photons are produced, within few tens fem- 
toseconds, at the same time. The probability of a spontaneous 
decay into a pair of correlated photons is usually very low, of 
the order of 1CP 9 or lower; therefore with typical pump power 
of the order of some milliwatts, the fluorescence emission lies 
at the levels of photon counting regime. Since the photons 
are produced in pairs and because of the energy and momen- 
tum conservation restrictions, the detection of one photon in 
a certain direction and with a given energy indicates the ex- 
istence of the pair correlated one, with definite energy in a 
well defined direction. This property allows to "heralding" 
the second photon of the pair once the first one is detected in 
a precise direction (and temporal and spectral window). This 
"heralded photon" was the state to be measured. In our set- 
up, after having eliminated the pump laser beam with a filter, 
the two photons of the pair are separated by means of a po- 
larizing beam splitter. The detection of one photon in one of 
the two conjugated directions was then used to open a win- 
dow of At = 4.9ns for detection in arm 2. This was realised 
by addressing the first detection signal as a start to a Time to 
Amplitude Converter, the signal from second detector (after a 
delay line) was then addressed to the same TAC as a stop and 
counted only if arriving in a window of At. 
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FIG. 3: Reconstruction of the photon distribution for the heralded 
single-photon state produced in spontaneous type II PDC. Inset: Ex- 
perimental frequencies /„ of no-click events as a function of the 
quantum efficiency r\ v for a PDC heralded photon state compared 
with the theoretical curve p v = 1 — f] v . 

The photodetection apparatuses were constituted by a sil- 
icon avalanche photodiode detector preceded by an iris and 
an interference filter (IF) at 702 nm, 4 nm FWHM, inserted 



4 



with the purpose of reducing the stray light. Both the detec- 
tors were silicon avalanche photodiode ones (SPCM-AQR-15, 
Perkin Elmer). The quantum efficiency of the "heralded pho- 
ton" detection apparatus (including IF and iris) was measured 
to be 20% by using the standard calibration scheme based 
on correlation properties of PDC emission (see fislo . Lower 
quantum efficiencies, needed for the reconstruction scheme, 
were simulated by inserting calibrated neutral density filters 
on the optical path. A comparison of the observed frequen- 
cies /„ with the theoretical curve (1 — rj v ) is presented in the 
inset of Fig. [5] The photon distribution has been reconstructed 
using K = 34 different values of the quantum efficiency from 
r\ v ~ 0.01% to r\ v ~ 20% with n v = 10 6 runs for each r\ v . 
Results at iteration i = 10 6 are shown in Fig. [3] As expected 
the PDC heralded photon state largely agrees with a single 
photon Fock state. However, also a small two photons com- 
ponent and a vacuum one are observed. The pi contribution is 
expected, by estimating the probability that a second photon 
randomly enters the detection window, to be 1.85% of pi, in 
agreement with what observed. A non zero po is also expected 
due to background. This quantity can be evaluated by mea- 
suring the counts when the polarization of the pump beam is 
rotated by a A/ 2 wave plate to avoid generation of parametric 
fluorescence. In this case as well, our estimate, (2.7 ± 0.2)%, 
is in good agreement with the reconstructed po- 

As a second example we have reconstructed the statistic of 
a strongly attenuated coherent state, which has been produced 
by a He-Ne laser beam attenuated to photon-counting regime 
by insertion of neutral filters. The same silicon avalanche pho- 
todiode detector of the previous case was used here as well. 
The counts were measured in about 400 ns window obtained 
by gating the photo-detector with a periodic signal (10 kHz 
rate). It must be noticed that in this case we do not have inter- 
ference filters or irises in front of the detector and all the other 
attenuations can be included in the generation of the state (i.e. 
they contribute to the absorbtion together with neutral filters): 
thus the highest quantum efficiency is assumed to be 66% as 
declared by the manufacturer data-sheet for the photodetector. 
The reconstructed distribution, with K = 15 different values 
of the quantum efficiency from r\ v ~ 0.1% to r\ v — 66% 
with n v = 10 6 runs for each r\ v , agrees well with what ex- 
pected for a coherent state with average number of photons 
\a\ 2 ~ 0.02. In Fig. 0]are shown both the frequencies /„ 
as a function of r\ v compared with the theoretical prediction 
p u — exp{— 77^ I ex 1 2 } ~ 1 — 77^ I « 1 2 and the reconstructed pho- 
ton statistic. Finally, we briefly acknowledge that a compa- 
rable result was also obtained with a very strongly attenuated 
thermal state, i.e. light emitted by a thermal source, a tungsten 
lamp attenuated by neutral density filters. 



IV. EXPERIMENTAL DATA IN PULSED REGIME 

In the pulsed domain, we have measured three different op- 
tical states generated starting from the third harmonics (349 
nm, 4.45 ps) of a cw mode-locked Nd:YLF laser regenera- 
tively amplified at a repetition rate of 500 Hz (High Q Laser). 
The general experimental setup is sketched in Fig. [5] For all 




FIG. 4: Reconstruction of the photon distribution for a weak coher- 
ent state. In the inset the experimental frequencies /„ of no-click 
events as a function of the quantum efficiency r\ v for a PDC heralded 
photon state are compared with a typical curve for a weak coherent 
beamp„ ~ 1 — r/ v \a\ 2 . 



the measurements, the light was delivered to a photo multi- 
plier tube (PMT, Burle 8850) through a multi mode fiber (100 
pm core diameter). Although the PMT has the capability of 
counting the number of photoelectrons produced by one or 
more photons 0], for the present application we used it in a 
Geiger-like configuration, by setting a threshold to discrimi- 
nate on/off events. Furthermore we take advantage of the lin- 
earity of the mean anodic charge, A c , see Inset in Fig. as 
a function of the mean energy of the measured light 1 20] to 
obtain the values of r) u . In fact, if we set A c — for 77 = 
and A c — (j4 c ) max for r] = rjp (nominal quantum efficiency 
of the PMT), for all the intermediate quantum efficiencies ob- 
tained by attenuating the light with neutral filters, we have 
V = ( 7 1p/ {A c )max)A c . This procedure allow us to vary con- 
tinuously the quantum efficiency. 
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FIG. 5: Setup for the generation of pulsed optical states. Starting 
with the third-harmonics of a Nd:YLF laser, we obtain: a) gaussian 
state; b) thermal state; c) multithermal state. D, rotating ground-glass 
diffuser; BBO, type I nonlinear crystal; F, variable filter; MF, mul- 
timode fiber; PMT, photomultiplier detector; GI, gated integrator; 
MCA, multichannel analyzer. 

The first measurement was performed on the pulse emerg- 
ing from the laser source. Due to the pulsed nature of the 
source, we do not expect to recover a true Poissonian statis- 
tics as in the cw measurement described above. Rather, we 
expect a Gaussian distribution of the form |21] 



1 



Qn,G 



■. exp 



(n — N) 2 
' 2{N + <t 2 



(9) 



which takes into account the presence of noise; N is the pho- 
ton mean value and a 2 /N can be taken as a measure of the 
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deviation from Poissonian statistics. In Fig. 0a) we show the 
photon distribution g n , reconstructed at the i = 50000 itera- 
tion of the ML algorithm, along with the best fit obtained with 
the model (|9j (fitting parameters N = 4.88 and a 2 = 0.63). 
The inset of the figure compares the experimental frequency 
/„ data (K = 37 values of rj, n v = 10 4 runs for each rj) as a 
function of r\ v with the theoretical values calculated through 
and the parameters given by the fit of the photon distribu- 
tion. Both the reconstructed distribution and the experimental 
frequencies agrees very well with the above Gaussian model. 
The fidelity of the reconstruction is G ~ 0.998. Using the 
estimated value of a 1 a deviation of about 13% of the laser 
photon number distribution from the Poissonian statistics can 
be derived. 
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FIG. 6: Histogram of the detected anodic charge of the PMT (data 
corresponding to the measurement of a gaussian field having N = 
4.88 with rj = 0.20). The arrow indicates the threshold value used to 
discriminate on/off events. Inset: linearity test for the PMT detector: 
mean anodic current as a function of the transmittance of calibrated 
filters. 



where /i is the number of temporal modes. The photon dis- 
tribution reconstructed at the i — 1500 iteration, is shown 
in Fig. c ) along with the best fit of the data using (^J 
(JV = 6.17 and /j, — 5); the fidelity of the reconstruction is 
given by G ~ 1. In the inset of the figure we show the ex- 
perimental frequency f v data (K = 18 values of rj, n u = 10 
runs for each rj) and their theoretical values as calculated ac- 
cording to Q. As a comment to the experimental results in 
the pulsed regime, we note that the best reconstruction of the 
photon distribution is achieved at a different number of itera- 
tions for the three different measured optical states, and that 
the absolute error e does not approach the same value. This 
is due to the presence of excess noise in our measurements, 
since the stability and the repetition rate of our source (500 
Hz) limits to n u ~ 10 4 the number of runs for each value 
of the quantum efficiency 1 13]. The choice of the best itera- 
tion to stop the algorithm is driven by the possibility to fit the 
distribution with a suitable model. We stress that there was 
no a-priori decision in choosing a Gaussian distribution for 
case a) or of a multithermal distribution for case c), but, on 
the contrary, we followed the a-posteriori observation that no 
other distribution could fit equally well the reconstructed data. 
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A second measurement was performed on the laser pulse 
diffused by a moving ground glass. When the photons are 
collected from within an area of spatial coherence, the sys- 
tem acts as a pseudo-thermal source, whose photon number 
distribution satisfies 

N n 

en,T - {N + l)n+ i ■ ( 10 ) 

Figure0b) shows the photon distribution g n , as reconstructed 
at the i = 400 iteration and the best fit of the data with (II Oi 
(JV = 5.33); the fidelity is given by G ~ 0.995. The inset of 
the figure contains the experimental frequency f u data (K — 
24 values of 77, n v — 10 4 runs for each rj) and their theoretical 
values as calculated from 0. 

The last measurement was performed on the blue portion 
(420 nm) of the down conversion fluorescence produced by a 
type I BBO crystal (10 mm depth, cut at 34 deg) pumped by 
the laser pulse. The pump, incident orthogonally to the crystal 
face, had an intensity ~ 60 GW/cm 2 . In this experimental 
condition we expect a coherence time of the generated field 
of ~ 1 ps, that corresponds to measuring a convolution of 
4-5 temporal modes 1 19]. The photon number distribution is 
expected to be a "multithermal" distribution of the form 

(n + n-l)\ 




FIG. 7: Reconstructed photon distribution (black bars) and best the- 
oretical fit (grey bars) for three different states in the pulsed regime: 
a) laser pulse, b) diffused laser pulse, c) multimode state produced 
in type I PDC. Insets: Experimental frequency /„ data in function of 
the quantum efficiency r\ v and theoretical model for each one of the 
states. 



Qn,M 



n\(n - 1)!(1 + JV/»"(1 + ri/NY 



(11) 
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V. CONCLUSIONS 

In summary, we have presented experimental results on 
the reconstruction of the photon distribution based on on/off 
detection at different quantum efficiency followed by a 
maximum-likelihood iterative algorithm based on the theoret- 
ical analysis presented in 11311 . 

Our results concern single-photon (PDC heralded) and 
weak coherent states in the cw regime, as well as to coher- 
ent, thermal and multi-thermal states in the pulsed regime. 
The stability and the good accuracy shown in the reconstruc- 
tion of these states, together with the simplicity of the method, 
clearly demonstrate the interesting potentialities of this tech- 
nique, suggesting relevant future applications ranging from 
studies on quantum optics, foundations of quantum mechan- 



ics, quantum information and quantum metrology. Some ap- 
plications in these fields are now under realisation in our lab- 
oratories (reconstruction of further quantum optical states, of 
entangled states 1 22], characterization of a high spectral selec- 
tion heralded photon source 12311 . etc.). 
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